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ABSTRACT 

We report the first counts of faint submillimetre galaxies (SMG) in the 870-/j,m band 
derived from arcsecond resolution observations with the Atacama Large Millimeter 
Array (ALMA). We have used ALMA to map a sample of 122 870-^m-selected sub- 
millimetre sources drawn from the 0.5° x 0.5° LABOCA Extended Chandra Deep 
Field South Submillimetre Survey (LESS). These ALMA maps have an average depth 
of <787o^m ~ 0.4 mJy, some ~ 3x deeper than the original LABOCA survey and 
critically the angular resolution is more than an order of magnitude higher, FWHM 
of ~1.5" compared to ~ 19" for the LABOCA discovery map. This combination of 
sensitivity and resolution allows us to precisely pin-point the SMGs contributing to 
the submillimetre sources from the LABOCA map, free from the effects of confusion. 
We show that our ALMA-derived SMG counts broadly agree with the submillimetre 
source counts from previous, lower-resolution single-dish surveys, demonstrating that 
the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. 
The difficulty which well-constrained theoretical models have in reproducing the high- 
surface densities of SMGs, thus remains. However, our observations do show that all 
of the very brightest sources in the LESS sample, S , 870/jm ^ 12mJy, comprise emission 
from multiple, fainter SMGs, each with 870-^m fluxes of ^ 9 mJy . This implies a nat- 
ural limit to the star-formation rate in SMGs of < 10 3 M Q yr^ 1 , which in turn suggests 
that the space densities of z > 1 galaxies with gas masses in excess of ~ 5 x 10 10 M Q is 
£ 10" 8 Mpc- 3 . We also discuss the influence of this blending on the identification and 
characterisation of the SMG counterparts to these bright submillimetre sources and 
suggest that it may be responsible for previous claims that they lie at higher redshifts 
than fainter SMGs. 
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1 INTRODUCTION 

The fir st deep surveys for extragalactic submillimetre 
sources (|Smail et al.lll997l :l Hughes et a l, 1998; Barger et al.l 
Il998h uncovered high number densities of submillimetre 
sources at mjy-flux limits and subsequent spectroscopy 
determined a median redshift of 2 ~ 2.5 for the radio- 
detecte d subset of the subm illimetre galaxy (SMG) popu- 
lation (|Chapman et al.ll2005h . At these high redshifts, the 
submillimetre fluxes of these sources correspond to far- 
infrared luminosities of >10 12_13 Lq, placing them in the 
ultra-luminous or hyper-luminous infrared galaxy (ULIRG, 
HLIRG) classes. These joint 850-/im and radio-selected sam- 
ples remain the best-studied SMGs and it has been claimed 
that they host up to half of the star formation occur- 
ring at z ;>2 (e.g.. |Hughes et all 1 19981 ; iBlain et all 1 19991 ; 



IChapman et al. 12005 ) and may be linked to QSO activity 
and the formation of massive galaxies a t high redshift (e.g., 
Swinbank et~ai1 120061 ; lAlexander et all 120081 ; iHickox et ail 



2012t ). If true then SMGs are an essential element in models 



of galaxy formation. In fact the first theoretical attempts 
to reproduce basic properties of SMGs, in particular the 
850- fim number counts, required radical alteration of the 
prescription for st arbursts in well-con s trained galaxy forma- 
tion models (e.g.. iBaugh et al.l I200I iGranato et al.| [2006), 
demonstrating the potential power of SMGs as a constraint 
on galaxy evolution theories. 

One concern about the use of the 850-/im number 
counts as a fundamental constraint on galaxy formation 
models is that these are derived from low spatial resolution 
(typically - 15-20" full width at half maximum, FWHM), 
single-dish surveys. This low resolution means that it is 
possible that several faint sources within a beam will ap- 
pear as a single brighter source, changing the shape of 
the number counts, most critically by potentially produc- 
ing a false tail of bright sources. A number of attempts 
have therefore been made to obtain high angular resolution 
continuum imaging through interferome tric observations of 
indiv i dual submillimetr e sources (e.g., iDannerbauer et al.l 
|2002|; IWang et aill201ll) and (nearly) flux-limited samples 



d Younger et al.ll2007l . l2009l ; ISmolcic et aT]|2012l ; iBarger et all 
2012). These observations have indeed shown that a number 
of bright submillimetre sources actually comprise emission 
from multiple SMGs. However, the conclusions from many 
of these studies have been weakened by a number of factors. 
In particular, both the modest numbers of sources studied 
and the fact that the discovery surveys underlying these 
studies are typically shallow or restricted to small areas, 
has meant it has not been possible to conclusively test the 
shape of the bright-end of the submillimetre source counts. 
Secondly, many of the follow-up observations of 870 /im- 
selected submillimetre sources have been carried out at 
longer wavelengths (typically > 1.2 mm) . This has led to am- 
biguous results, especially when comparing the single-dish 
and interferometer-based fluxes for sources, and so limits 
the conclusions that can be drawn about their multiplicity. 
However, on e particularly notew orthy stu dy is that recently 
published bv lBarger et al.l l|2012l ) (see also lWang et al.ll201ll ) 
which used the Submillimeter Array (SMA) at 850 )im 
to observe sixteen 850 ^tm-selected submillimetre sources 
from a SCUBA-2 survey of the Gr eat Observatories Ori- 
gins Deep Survey North (GOODS-N; lGiavalisco et alj|2004 



I Wang et al.ll20o3 ). This wavelength-matched study showed 
the best evidence yet for an increased incidence of multiple 
SMGs in submillimetre sources at bright 850-/im fluxes. 

The issue of reliable SMG counts, needed to robustly 
constrain the theoretical models, is obviously an area where 
the Atacama Large Millimeter Array (ALMA) will have sig- 
nificant impact. In 2004 we therefore started planning a sur- 
vey to provide a large, flux-limited sample of submillimetre 
sources over a wide area in a field with excellent visibility 
from ALMA. The field chosen was the 0.5° x 0.5° Extended 
Chandra Deep Field South (ECDFS) which has the most 
extensive multi-wavelength coverage of any large-area ex- 
tragalactic region in the southern hemisphere. The result 
was th e LABOCA ECDFS submillimetre survey I Weifi et all 
(2009) (LESS and W09 hereafter), which obtained a deep, 
o"s70;jm ~1.2mJy, homogeneous 870-/xm map of the full 
ECDFS detecting 126 submillimetre sources. 

As the next step, in ALMA Cycle we observed 122 
of the 126 submillimetre sources from LESS using ALMA 
in its compact configuration. These data and the result- 
ing SMG catalogue will be presented in Hodge et al. (2012, 
in prep). Critically, these observations were carried out at 
the same wavelength as the LABOCA survey. The resulting 
ALMA maps yield unambiguous identifications for a large 
fraction of the submillimetre sources, directly pin-pointing 
the SMG(s) responsible for the 87 0- /xm emission to within 
<0.3"; (e.g. ISwinbank et al.ll2012T )). The spatial resolution 
achieved by our observations (~ 1.5" FWHM) corresponds 
to an order of magnitude improvement over the single-dish 
LABOCA survey. It thus provides an ideal data set to deter- 
mine the influence of multiplicity on the form of the 870-/im 
SMG counts, as require d to enable a reliab le comparison to 
model predictions (e.g.. IBaugh et a"lll2005l ). 

In this paper we analyze these ALMA maps to derive 
number counts for SMGs and compare these to both previ- 
ous source counts from single-dish submillimetre surveys and 
to predictions from theoretical models. We adopt a cosmol- 
ogy with S1 A =0.73, Q M =0.27, and H =72kms _1 Mpc" 1 
in which a scale of 1" corresponds to a physical separation 
of ~ 8.4 kpc at a redshift of 2 = 2. 



2 OBSERVATIONS AND CALIBRATION 

Of the 126 submillimetre sources detected in the LESS sur- 
vey, 122 were observed during Cycle using the band 7 re- 
ceivers in ALMA's compact array configuration. The avail- 
able 8-GHz bandwidth was centered at an observed fre- 
quency of 344 GHz (i.e. 870 fim) and we employed a dual po- 
larization setup. This choice of observing frequency enables 
a direct comparison to the flux densities of the submillimetre 
sources measured in the original LABOCA survey. The ob- 
serving campaign was carried out in eight observing blocks 
(ALMA measurement sets; MS hereafter) between 2011 Oc- 
tober 18 and 2011 November 3. Typically 15 antennas were 
available for each block. 

The quasar B 0402-362 (J0403-360) was used for 
phase calibration and Mars, as well as Uranus, have been 
used to calibrate the absolute flux scale. Bandpass cal- 
ibration was generally performed using observations of 
B 0537-441 (J0538-440). Each science field, centered on 
the catalogued position of a given LESS source from W09, 
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Figure 1. Examples of the 870-/^m ALMA continuum maps towards eight of the submillimetre sources from the LESS survey. In each 
map we identify all of the sources with S / N > 3.5 it (squares). The ALMA data unambiguously locates the SMGs to a precision of < 0.3" 
and to flux limits of ~2mjy (~ 3.5-7 a). The upper row shows the maps of the brightest LESS sources in our main sample (see text for 
details) while the lower row shows maps that are randomly drawn from each quartile of the LESS flux distribution. Positive (negative) 
contours on each map are shown in black (white) and start at (— )3<r and are incremented by (— )2a. The 1-c noise in the map is shown 
in the bottom left corner of each panel. Each map is 25.6" across and we show the primary beam (dotted circle) encompassing the radius 
at which the ALMA antenna sensitivity drops to 50%. 



was observed for a total of ~ 120 seconds. The data were pro- 
cessed with the Common As tronomy Software Application 
(CASA; iMcMullin et al.ll2007f ) and imaged using the CLEAN 
algorithm within CASA. A detailed description of the raw 
data and its calibration as well as imaging will be presented 
in an upcoming publication (Hodge et al. in prep.). 

The field of view - defined as the FWHM of the ALMA 
antenna reception pattern around the phase center and re- 
ferred to as primary beam in the following - is 17.3" in 
diameter. Q Each map has a pixel size of 0.2" and a total 
extent of 128 pixels in each dimension, sufficient to cover the 
primary beam and encompass the error-circles of the sub- 
millimetre sources from the LESS maps, ^ 5" (W09), even in 
confused situations. The average root mean square (rms, a) 
of the background noise in the maps is a ~ 0.4 mjy - a fac- 
tor ~ 3 deeper than the original LABOCA observation. Us- 
ing natural weighting we achieve a typical restoring (clean) 
beam of ~ 1.8" x 1.2", although a small number of low el- 
evation (<C30deg) observations lead to much larger beam 
ellipticities and the corresponding image products are typi- 
cally much noisier than our median maps, producing a tail in 
the noise distribution extending beyond 0.6 mjy. In the fol- 

Accordingly, the flux density at a given position in the resulting 
map can be corrected by multiplication with the corresponding 
factor derived from an inverse 17.3" FWHM Gaussian. We will 
note in the following when primary-beam corrected fluxes are 
used. Note that the FWHM used, based on actual beam measure- 
ments is slightly smaller than theoretically expected for a 22-m 
antenna (ALMA ARC, private communication). 



lowing - unless explicitly stated otherwise - we will focus on 
the sub-set of 88 "best" maps selected from two important, 
but not mutually exclusive, selection criteria: beam-axial ra- 
tio < 2 and rms noise level a < 0.6 mjy0 The distribution of 
targets between and within MS was chosen so that problems 
with the observations of any particular MS would not bias 
our sample and hence our "best" sample represents a ran- 
dom sampling of the LESS catalogue, yielding an unbiased 
view of the properties of submillimetre sources as a function 
of 870-^im flux. Figure Q] shows examples of the calibrated 
and cleaned maps used in this study. 

3 IDENTIFICATION OF ALMA SOURCES 

3.1 Source extraction and characterization 

In order to detect SMGs in our calibrated and cleaned maps 
we use an automated scheme (described in detail in Hodge 
et al. in prep.). Our iDL-implemented source extraction soft- 
ware first identifies individual signal peaks above a 2.5-ct 
threshold which are used as the basis to model the emission 
in that region using a Metropolis-Hastings Markov chain 

2 Note that we analyze those maps that do not comply with our 
selection criteria in the same way as the 88 "best" maps, including 
the identification of sources as described in the following. Bright 
sources detected at sufficient significance in those maps form part 
of a supplementary ALMA catalogue which is not used in this 
paper unless explicitly stated. 
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Monte Carlo (MH-MCMC) algorithm to determine the best 
six parameter fit for an elliptical Gaussiarjf] to describe the 
underlying flux distribution within a 2" x 2" region. This 
is large enough to recover any extended sources but small 
enough to resolve double sources. Whilst we attempt a full 
six parameter fit, at the resolution of our survey (~ f .5" 
corresponding to physical scales of ~ 12 kpc at z > 1), many 
sources are unlikely to be resolved. We therefore repeat the 
fitting process using a simple elliptical point-source model 
with only three free parameters and source extents as well 
as orientation fixed to the synthesized beam parameters. We 
find that the peak flux densities from the simple fit generally 
agree with the integrated flux densities given their full fit- 
ting errors, we find a > 2cr integrated flux e xcess in just two 
source s, one of which has been discussed bv lSwinbank et aU 
(2012). For all other sources we therefore adopt the peak 
flux density from the three parameter model as our best 
estimate of the source flux. 

For each parameter in the fit its posterior distribu- 
tion determines the fitting error, after correcting for auto- 
correlation in the Markov chain. These uncertainties are 
taken into account when determining the full measurement 
error for the integrated source flux density. The recipes 
we follow to determine the full uncertainty for elliptical 
source fits in the presence of correlate d nois e in radio 
maps have been motivated by Condon! (11997 ) (see also 
Windhorst et al.lll984l;lHopkins et al.ll2003l . lSchinnerer et all 
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Figure 2. Results from artificial source simulations to test the re- 
liability of our source extraction procedure. After removing > 2.5a 
sources identified by our automatic source finder we insert in each 
map five well-separated sources at random positions within the 
primary beam, then run our source extraction code and repeat 
this process 16 times. Red squares denote the fraction of inserted 
sources that are recovered within a given narrow signal-to-noise 
ratio (SNR) bin while blue circles indicate the fraction of ex- 
tracted sources which are unassociated within 0.8" with an in- 
serted source. Above the 3.5-cr threshold (vertical dashed black 
line) - adopted as the reliability limit for our analysis - our cata- 
logue is ~ 99% complete and has a false detection rate of ~ 1.6%. 



3.2 Source extraction efficiency and flux recovery 

Our goal is to construct a catalogue which is deep, i.e. in- 
cludes as many real faint sources as possible, but has a very 
low spurious source fraction. To determine the search param- 
eters necessary to achieve this we use a suite of simulated 
ALMA maps. First we prepare cleaned map which are our 
actual ALMA maps with all sources > 2.5a identified by our 
automatic source finder removed. Within the primary beam 
area of each cleaned map we insert five sources separated 
from each other by at least two synthesized beam widths. 
The inserted sources cover a wide range in significance (~ 2- 
20 a) and follow a steeply declining flux density distribution. 
We then run our code and extract all sources in the simu- 
lated maps, in exactly the same fashion as the science maps. 
We then determine the recovery rate of our inputs sources 
and how many spurious sources, not associated with any of 
our model sourcefl are found. This process is repeated 16 
times per map. Thus, in total we insert 7,400 sources into 
the 88 ALMA maps. 

Figure [2] shows the results of our simulations for the 
fraction of recovered sources and the fraction of detected 
sources which are spurious, both as a function of mea- 
sured peak signal-to-noise ratio. The cumulative distribu- 
tion shows that the source extraction recovers ~ 99% of all 

3 These six parameters are peak flux density, pixel position of 
the peak, minor axis extent, axial ratio and position angle 

4 Below our initial CLEAN threshold our simulated sources are in- 
serted after convolution with the dirty beam instead of the CLEAN 
beam used for brighter sources. The dirty beam side lobes can 
then be boosted by noise peaks to appear as spurious sources 
above our detection limit. 



> 3.5a sources, while a sample selected above this thresh- 
old contains only ~ 1.6% spurious detections. We therefore 
adopt a 3.5-cr detection limit for the main (i.e. most robust) 
sample used in our analysis. We expect that this sample will 
contain less than two spurious sources and will fail to include 
only one intrinsically > 3.5-cr source. We also employ the re- 
covered and spurious fractions as a function of detection 
significance to correct our measured counts in our analysis 
below. These corrections become significant only below the 
3.5-cr catalogue limit and we highlight the flux regime most 
affected by those corrections when discussing the counts. 

As a final test we compare the total fluxes of the SMGs 
detected by ALMA to the deboosted fluxes measured for the 
submillimetre sources in W09. To achieve this we sum the 
primary-beam-corrected ALMA fluxes of all sources above 
a given detection significance, within the ALMA primary 
beam area, weighted by the LABOCA beam to calculate 
the total flux that would have been seen by LABOCA at the 
submillimetre source positions from W09. To our 3.5-cr sig- 
nificance threshold the resulting median ALMA/LABOCA 
flux density ratio and bootstrap error is 0.83±°;°*. K the 
LABOCA flux scale and the flux deboosting of the LESS 
sources are accurate, this suggests that a contribution from 
additional fainter SMGs is required to recover the total flux 
density within the LABOCA beam. Integrating the flux in 
sources down to a 3-cr significance limit results in a median 
flux ratio of 1.05 ±0.09, consistent with unity. Given the in- 
creasing number of spurious detections at lower detection 
significance we will, however, retain a 3.5-cr limit for our 
main sample and discuss the comparison of flux scales in 
more detail in Hodge et al. (in prep.). 

In total we detect 99 individual SMGs at > 3.5 a in the 
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88 maps in our main sample. Of these maps, 69 show at 
least one SMG, 19 maps exhibit two SMGs and four maps 
have three SMGs within the primary beam. Hence 19 ALMA 
maps (~22%) do not contain a >3.5ct source. The asso- 
ciated LABOCA submillimetre sources have a median de- 
boosted flux of 4.5 mjy and a median detection significance 
of 3.75 a, making the sources amongst the faintest in the 
LESS survey. The resulting full ALMA SMG catalogue and 
the exploitation of this catalogue and maps to investigate the 
multi-wavelength properties of the SMGs will be presented 
in upcoming publications (Hodge et al. in prep.; Simpson et 
al. in prep.). 



Table 1. Parameterized fits to the differential counts of 870- /^m 
ALMA-detected SMGs 



dN/dS v N* S* a 

[1/N*] [mJy-Mcg- 2 ] [mjy] 



a) "e-'Ws; 27 0.9 -3.6 

b) (fr)~°\ S v < St 

(§£■) S ,S v ^S* 1564 8 2.1 7.8 



4 RESULTS AND DISCUSSION 

The ALMA SMGs we have identified can be used to estimate 
the 870-/im source counts free from the influence of blending. 
In the following we describe the derivation of these counts 
and compare the results to those from previous single-dish 
surveys. 



4.1 Derivation of ALMA 870 (im source counts 

We are interested in the differential number counts of SMGs 
as a function of flux. Our ALMA survey covers submillimetre 
sources above the flux limit of the LESS discovery survey 
catalogue (W09) so that the effective ALMA survey area, 
Aalma, above this selection limit is given by 



Aalma 



Amaps (ALMA) 

(W09) 



x Aless, 



(1) 



so urc cs 



where Af maps (ALMA) =88 is the number of ALMA maps 
used in this study, Aless — 0.35 deg 2 is the LESS survey 
area and iV SO u rccs (W09) = 126 is the total number of sub- 
millimetre sources in the LESS catalogue. 

Necessarily, this area is only correct in the flux regime 
covered by the LESS surveyQ Below the LESS flux limit the 
source counts derived from our ALMA maps must be con- 
sidered biased since the observations were typically taken in 
the vicinity of a brighter submillimetre source and so are not 
necessarily representative of the fainter source population. 
This restriction already applies to the faint-end of our main 
sample of > 3.5-cr ALMA SMGs (which are below the LESS 
catalogue flux limit) and we neither attempt any interpre- 
tation of these values nor apply further assumptions. For 
completeness and purely informational value, we addition- 
ally derive differential counts for all > 2.5-ct ALMA SMGs. 

Given our findings in £13.21 the source sample and hence 
the counts are affected by our signal-to-noise ratio (SNR) 



5 Strictly speaking, this area is only valid for sources in the phase 
center of a given ALMA map since the sensitivity monotonically 
drops to 50% at the primary beam radius. Nevertheless, for the 
interesting regime of our analysis - above the LESS survey flux 
limit — even at half the phase center sensitivity sources would be 
detected at a > 3.5 a level and hence included in our main sample 
even if residing right at the edge of our field of view. 



Notes: Best-fit parameters to describe our data: a) Schechter func- 
tion, b) broken power law. Both fits are applied only to the data 
above the LESS catalogue limit and provide a good representation 
of the steepness of the bright end (> 8 mjy) of our source counts. 
We urge caution when extrapolating these functions, particularly 
the Schechter parameterization, to much fainter fluxes. 

selection and so we use the recovery and spurious frac- 
tion rates derived from our simulations to correct our ob- 
served counts. However, we stress that these extraction bi- 
ases barely affect our main sample and even less so the flux 
regime above the LESS detection limit. Our paramctcri- 
zations (from Figure [2)| of the fractions of spurious detec- 
tions, /spurious (SNR), and sources recovered in our simula- 
tions, /recovered (SNR), provide us with the probability that 
an SMG with a given SNR is spurious and also the like- 
lihood that SMGs of the same SNR are missed. For each 
SMG we therefore derive an individual source probability 



ofp(SNR) = 2-f r< 



i(SNR)-t, 



9 (SNR). Corrected 



differential source counts are then given by the total of all 
source probabilities in a given flux bin and normalized using 
Equation ([T}. 

Figure [3] shows the resulting differential as well as the 
cumulative number counts along with the corresponding cor- 
rected values. The differential source counts derived for our 
main sample (> 3.5 a) are, as expected, only affected by the 
bias-corrections at the lowest flux densities while the cumu- 
lative counts are even less affected. For flux density bins con- 
taining < 3.5-cr sources - far below the LESS survey limit - 
we show only bias-c orrected values. We also show the best-fit 
Schechter function l|Schechterll 19761 ) and a broken power-law 
to describe our differential counts above the LESS catalogue 
limit and summarize all these fit parameters in Table [1] 



4.2 The absence of the bright SMG population 

The most surprising result of our source counts is a clear 
break at the bright end caused by a lack of bright SMGs in 
our ALMA maps. None of our maps detect an SMG with 
a flux ^9mJy despite twelve sources with 870- ^tm fluxes 
> 9 mjy in the LESS survey. As Figure Q] shows, in at least 
four of these cases, the submillimetre sources comprise mul- 
tiple (fainter) SMGs. This is particularly clear in the very 
brightest LESS sources (^12mJy) where we detect mul- 
tiple high-significance (^> 6a) SMGs in each map. In the 
remaining cases, a single ALMA SMG is detected above 
3.5 a, although in all cases, the flux density of this SMG 
significantly under-predicts the LABOCA flux. This short- 
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Figure 3. Left: The primary-beam corrected differential 870-fim source counts from our ALMA survey. Above the LESS survey limit 
(depicted with a shaded margin representing the LABOCA rms) the counts of ALMA sources should be a true representation of the 
underlying SMG population. We also plot counts below this limit for our robust ALMA sample (all bins depicted by filled circles comprise 
> 3.5-cr detections) which is largely uncontaminated by spurious sources while showing a high detection efficiency (see Fig.fJJ. In addition 
we extend the counts to > 2.5 cr significance applying the larger corrections necessary for the fractions of spurious and undetected sources. 
Fainter than ~ 9 m Jy our data are slightly lower but in reasonable agreement - within the error margins - with the scaled fit (dot-dashed 
line) to the number counts from the SHADES survey llCoppin et al.l l2006'). However, brighter than ~ 9 mjy we find a steeper decline 
than shown in the single-dish counts (a Schechter parameterization of the differential counts resulting from a P(D) analysis of the LESS 
map by W09 similarly over-predicts the bright counts). Submillimetre sources brighter than this limit were detected by the LESS survey 
(deboosted values; blue filled squares) but no comparably bright SMG is observed by ALMA. Instead the ALMA maps of the brightest 
LESS sources show an increase in source multiplicity. To parameterise our counts we fit two parametric models to the data above the 
LESS survey limit and overplot these here: A Schechter function and a broken power-law (see Table [TJ shown by red dashed and solid 
line s respectively . Rig ht: Corresponding cumulative counts of > 3.5-cr ALMA sources compared to the deboosted LESS results and 
the iBaugh et al] l|2005l ) model predictions. The latter are constrained by observed data from single-dish submillimetre surveys and the 
predicted fraction of bursts is shown in addition to the total counts. Scaled to the LESS data the model is comparable to our results 
above the selection limit while the clear mismatch at the bright end due to the absence of corresponding ALMA SMGs is apparent. The 
cumulative counts are corrected for the fraction of missing/spurious detections, but this has a minor effect to our main sample due to 
its high cumulative detection efficiency. 



fall in flux could arise either from our resolving-out extended 
emission or from the presence of several ~ 1-1.5 mjy SMGs, 
which lie below our detection threshold. Nevertheless, an im- 
portant result of our survey is that the brightest submillime- 
tre sources (> 9 mjy) in single-dish surveys likely comprise 
multiple, fainter SMGs. 



SMGs with 870-/im fl uxes of > 9 m jy are 
likely to be HLIRGs (d Rowan-Robinsonl l200d : 
iRowan-Robinson fc Wanel l2010h ) if they lie at z > 1. 
The lack of large numbers of such bright sources in our 
sample then implies a natural limit to th e star-formation 
rate in an SMG of £ 10 3 M Q yr" 1 (for a ISalpeterl (|l955l ) 
IMF). This maximal star-formation rate (SFR) is driven 
by the ratio of the mass of the available gas reservoir an d 
the free-fall time of the system l|Lehnert fc Heckm an 1996). 
For a free-fall time of ~ 50 Myrs as found by Ke nnicuttl 
( 1998) for local starbursts, the implied gas mass limit 



is g 5 x 10 10 Mq. This is comparable to the l imiting gas 
mass for SMGs found by iBothwell et al.l (|2012T ). Thus the 
absence of very bright SMGs in the LESS survey area 
therefore suggests that high-redshift galaxies with gas 
masses significantly above 5 x 10 10 M© have space densities 
of ^ 10~ s Mpc~ 3 . For normal star forming z < 3 galaxies 
iKarim et al] l|201ll ) suggest that their inverse free-fall time 
constitutes a potential upper limit to their specific SFR. 
Assuming a typical stellar mass of ~ 10 11 M for our SMGs 



(e.g., ISwinbank et al] |2GT2| ) also their specific SFR is in 
agreement with this upper limit. 

As a consequence of the absence of bright SMGs, our 
ALMA source counts show a sharp drop at ^ 8 mjy (Fig- 
ure [3]l . Parameterising the differential counts with a double 
power-law broken at this flux value, we find a factor of > 4 
difference in the power law indices of the two parts (see 
Table [TJ. We caution that our data coverage of the regime 
brighter than the LESS survey limit is really too sparse for a 
four parameter model fit. However, a simple power-law fit to 
the breaking point - effectively representing a step function 
- leads to the same result. Similarly, a Schechter function ap- 
pears to be an equally good description of our data above the 
LESS survey threshold, although we urge caution to extrap- 
olate this parametrization into the low flux regime where 
our data cease to be representative. We hence conclude that 
the broken power-law constitutes the best representation of 
our data. 

Although previous interferometric surveys of submil- 
limetre sources cannot be considered complete as they com- 
prise a complex mix of follow-up observations from hetero- 
geneous surveys and small survey fields, it is instructive to 
co mpare our findings t o these previous results. For exam- 
ple iBareer et ID (I2012T) obtained deep integrations of four 
JilOmJy SMGs in GOODS-N with the SMA at 860 pm 
and showed that at least two break up into multipl e com - 
ponents across ~5". Earlier results by IWang et al] (|201ll ) 
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also suggested that potentially ~ 30% of > 5-mJy 850-/im 
sources could comprise such multiple sy stems, potentially 
rising to > 90% above ~8mJy (see a l so I Wang et alJ 120071 : 
lYounger et"afll200Sl : ICowie et al.ll2009l ; ISmolcic et alj|2012h . 

These results support our interpretation that the num- 
ber of bright submillimetre sources £ 9 mjy from single-dish 
surveys have been substantially over-estimated, producing 
artificially high SMG number counts at the brightest fluxes. 
To highlight this we show in Figure [3] the results from W09 
and those from the SCUBA HAlf D egree Extragalactic Sur- 
vey (SHADES; ICoppin et ahlboOrj ) scaled to the LESS data 
as described by W09. We refer to both studies for an exten- 
sive comparison to other ~ 850-/^m single-dish surveys. It 
is noteworthy that even sophisticated methods to estimate 
the differential source counts from single-dish submillime- 
tre surveys - such as the probabilistic P{D) analysis pre- 
sented by W09 - does not remove the bright source excess 
(see Figure (3|. While a P(D) analysis, deriving the differen- 
tial counts solely from the flux distribution in the map and 
not based on individually extracted sources, would need to 
account for the realistic clustering statistics on small angu- 
lar scales in order to derive the true counts. Interestingly, 
our counts at intermediate flux levels, while not being sig- 
nificantly lower, are also not particularly higher than the 
single-dish results, indicating that below the break we also 
find fainter sources suffering from multiplicity effects. This 
suggests that those ALMA maps of fainter LABOCA sources 
which lack > 3.5-<r SMGs may be explained by the presence 
of multiple SMGs below our significance limit[f] Hence, the 
statistical prediction by W09 of only five LESS sources being 
spurious may still be valid. 

The high multiplicity of the brightest submillimetre 
sources also has wider implications. For example, it has been 
claimed that the brightest (most luminous) SMGs evolve 
more strongly than fainter systems and hence are preferen- 
tially found at the highest redshifts (e.g., Ivison et alJ 20021. 
20071 : IWall et all 120081 : iMarsden et al.ll201ll ; ISmolcic et all 
20121 ). This could now simply be explained by two effects: 
confusion, which nullifies the statistical techniques used to 
identify counterparts, meaning that the very brightest sub- 
millimetre sources lack obvious counterparts at other wave- 
lengths and as a result are associated with (undetectable) 
high-redshift sources; and the significant overestimation of 
the submillimetre fluxes for the identified counterparts, 
whose artificially enhanced radio/submillimetre flux ratios 
then mimic those expected for high-redshift SMGs. 

Similarly, the detection rate of 12 CO emission in the 
SMG survey of iBothwell et al l (|2012h declines with 850- /im 
flux. For example, from their sample of 40 SMGs, the eight 
of which are undetected (including 2/5 with S$so ^ lOmJy) 
in 12 CO have a median Ssso = 8.1 ± 0.7mJy, compared to 
Ssso =5.9±0.7mJy for the 12 CO detections (this is signifi- 
cant at 92% confidence level). This modest difference would 
be explained if a higher fraction of the brighter sources have 
their submillimetre fluxes boosted by emission from other 
sources projected along the line of sight withi n the beam. 
This scenario is consistent with the suggestion of lWang et al.l 



l|201ll ) that where several SMGs comprise a submillime- 
tre source, those individual components are not necessarily 
physically associated. 

The high multiplicity of bright submillimetre sources 
may also bias the form of the far-infrared-radio correla- 
tion, which is widely used to infer star formation rates 
out to high redshifts. For example, where the submillime- 
tre flux of a source is derived from single-dish photome- 
try, and in fact represents contributions from several SMGs, 
with the radio emission coming from a single source, then 
this will produce a systematic offset and scatter in the de- 
rived far-infrared-radio correlation. A number of studies 
e.g.Jlbar et al.ll2008l ; ISargent et aill2010l : llvison et al1l2010l : 



iBourne et al. 201 ll ) have attempted to trace the evolution of 
this relation over a range of redshifts, of whi ch several fo- 
cussed on submillimetre-b r ight sources (e .g., iKovacs et al.l 
120061 ; iMurphv et al] 120091 : iMurphvl I2009T ). The latter re- 
ported a radio excess for submillimetre-selected sources com- 
pared to local star-forming galaxies and - if true - this offset 
would increase further if only a single sub-component is as- 
sociated with the radio emitter^ However, we postpone a 
discussion of the radio properties of our ALMA SMGs to an 
upcoming publication. 

Finally, we also compare our SMG count s in Figure [3| 
to th e predicted counts from the model of iBaugh et al.l 
(2005) which, by design, match the single-dish observations. 
In comparison, consequently, our counts are significantly 
lower, particularly at the bright end. Since IBaugh et al.l 
(2005) could only reproduce the submillimetre counts by 
including a top-heavy IMF during bursts of star forma- 
tion, the tension created by our ALMA counts might aid 
the model to reproduce the true counts without resorting 
to a non-standard IMF. However, our data show no strong 
disagreement at intermediate fluxes so that the fundamen- 
tal problem for models of reproducin g the normalisa t ion o f 
the counts may still persist. Recently. lHavward et al.l (120121 ') 
predicted the submillimetre source counts by a combined 
semi-analytical/hydrodynamic approach without including 
a top-heavy IMF for bursts. Throughout the flux regime 
covered by our study, they predict a consistently high frac- 
tion (> 30%) of galaxy pairs, in agreement with our observa- 
tions. Still, their predicted bright SMG counts do not drop 
as sharply as our data indicates. 



5 SUMMARY 

We have presented source number counts derived from an 
870-/im ALMA survey of submillimetre sources from the 
870-/^m LABOCA survey of the Extended Chandra Deep 
Field South by W09. Compared to the parent survey, our 
ALMA maps are three times deeper and have an angular 
resolution which is an order of magnitude higher, allowing us 
to remove the influence of blending on the counts above the 
LABOCA detection limit. We find that our source counts 
are in broad agreement with those of the LABOCA sur- 
vey and previous literature results. However, brighter than 
~ 8 mjy our counts show a deficit of sources compared to 



6 Assuming standard dust properties a galaxy at 2 > 1 with a 
star-formation rate of 100 Mq yr — 1 would have an 870-/im flux of 
~ 1 mjy, below our detection limit. 



7 It should be noted that Bar ger et al. I j2012l) find that five of 
their SMA sources agree with the local relation. 
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those from single-dish surveys. This is caused by multiple 
SMGs, which are resolved at the ~ 1.5" resolution of our 
ALMA maps, being blended into single sources at the resolu- 
tion of the single-dish surveys. This trend has also been seen 
in recent studies of significantly smaller samples of submil- 
limetre sources. Our results suggest that multiplicity in sub- 
millimetre sources is significant at the brightest fluxes, but 
may also influence fainter submillimetre sources, ~ 4 mjy, 
from single-dish surveys. The absence of bright SMGs in our 
sample implies a limit to the maximum star-formation rate 
in an SMG of & 10 a M Q yr" 1 (for a Salpeter IMF), which 
in turn suggests that systems with gas masses in excess of 
~5x 10 10 M Q have space densities of <, 1CT 8 Mpc -3 . 
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